KEK-TH-1232, RBRC-756, YITP-08-74 



PROCEEDINGS 

OF SCIENCE 



Nucleon structure with dynamical (2+1)-flavor 
domain wall fermions lattice QCD 



Shigemi Ohta* 

Institute of Particle and Nuclear Studies, KEK, Tsukuba, Ibaraki 305-0801, Japan 
RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973 
Physics Department, Sokendai Grad. Univ. Adv. Studies, Hayama, Kanagawa 240-0193, Japan 
E-mail: shigemi . ohtaSkek . jp 

Takeshi Yamazaki^ 

Yukawa Institute for Theoretical Physics , Kyoto University, Kyoto 606-8502, Japan 
E-mail: tyamazakSyukawa . kyoto-u .ac.jp 

RBC and UKQCD Collaborations 

We report isovector form factors and low moments of isovector structure functions of nu- 
cleon from the coarse RIKEN-BNL-Columbia (RBC) and UKQCD joint dynamical (2+l)-flavor 
domain-wall fermions (DWF) ensembles. The lattice cut off is estimated at a^' = 1.7 GeV. The 
lattice volume is as large as 2.7 fm across. The dynamical strange mass is set at slightly heavier 
than physical, and degenerate up and down mass is varied correpsonding to pion mass of about 
0.67, 0.56, 0.42 and 0.33 GeV. We carefully optimize the nucleon source/sink separation in time 
to about 1.4 fm. Unexpectedly large finite-size effect in the axial charge is found. The effect 
scales with a single variable, the product nij^L of the pion mass nij^ and lattice spatial linear ex- 
tent L, and sets in at around nij^L = 5. We also discuss momentum-transfer dependence of the 
vector, induced tensor, axial-vector and induced pesudo-scalar form factors. In particular we ex- 
tract Dirac and Pauli mean-squared charge radii, magnetic moment, nNN coupling constant and 
pseudoscalar coupling gp of muon capture by nucleon. While the mean-squared charge radii 
are distant from the respective experiments, the other observables are in reasonable agreement 
with the experiments. From structure functions, fully non-perturbatively renormalized iso-vector 
quark momentum fraction, (x)„_j, helicity fraction, {x)\u-M^ and transversity, {y)su-Sch re- 
ported, as well as an unrenormalized twist-3 coefficient, d\ . The ratio of the momentum to helicity 
fractions, {x)„^ii/ {x)Au-Aih does not depend on light quark mass and agree well with the experi- 
ment. Their respective absolute values, fully renormalized, shows interesting trending toward the 
respective experimental values at the lightest light quark mass. 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker of the structure function part. 
^Speaker of the form factor part. 



© Copyright owned by the author{s) under the terms of the Creative Commons Attribiition-NonCommercial-ShureAHke Licence. 



http://pos.sissa.it/ 



(2+1 )-flavor dynamical dwf nucleon structure 



1. Introduction 

We report numerical lattice quantum chromodynamics (QCD) calculations of all the four 
isovector form factors and some low moments of the isovector structure functions of nucleon us- 
ing the lattice gauge ensembles [1] jointly generated by the RIKEN-BNL-Columbia (RBC) and 
UKQCD Collaborations with "2+1" flavors of dynamical domain-wall fermions (DWF) [2, 3, 4]. 

Four isovector form factors parameterize neutron j8 decay: the vector and induced tensor form 
factors from the vector current, 

{p\V+{x)\n) = Up [rnGv{q^)-qp.ax^GT{q^)] uj'^\ (1.1) 

where Gy is equivalent to Fi and Gt to F2/{2m!\i) in electromagnetic form factors, and the axial 
and induced pseudoscalar form factors from the axial current, 

{p\A+{x)\n) = Up [Y^r5GA{q^)+iq^r5Gp{q^)] Une'i'. (1.2) 

We use the Euclidean metric convention: thus q^ stands for Euclidean four-momentum squai^ed, 
and corresponds to the time-like momentum squared as = — < in Minkowski space. Here 
q = p„— pp is the momentum transfer between the proton (p) and neutron (n). In the limit \q\ -^0, 
the momentum transfer should be small because the mass difference of the neutron and proton 
is only about 1.3 MeV. This makes the limit q^ — > 0, where the vector and axial form factors 
dominate, a good approximation. Their values in this limit ai^e called the vector and axial charges 
of the nucleon: gy = Gv{q^ = 0) and gA = GA{q^ = 0)- As is well known, the vector charge 
is determined by the Cabibbo mixing angle dc'- gv = cos 0c- The axial charge receives strong 
interaction correction, and is very accurately known, gA = 1.2695(29) x gy [5] from neutron j8 
decay measurements. It is an interesting and important challenge for numerical lattice QCD to 
reproduce this value. 

Nucleon form factor experiments of course are not limited to neutron j3 decays. Vast data of 
lepton elastic scattering exist at wide range of momentum transfer values. On the lattice we can 
investigate nucleon form factors at relatively low momentum transfer such as ^ =~ In/Lx (0, 0, 0), 
(0,0,1), (0,1,1), (1,1,1), (0,0,2), characterized by the inverse of the linear lattice spatial extent L. 
From the dependence on momentum-trasfer we can extract mean-squared radii of the respective 
form factors, anomalous magnetic moment, nNN coupling, pseudoscalar effective coupling gp 
that enters the muon capture by nucleon, and so on. 

The structure functions are measured in deep inelastic scattering of leptons from nucleons, the 
cross section of which is factorized in terms of leptonic and hadronic tensors, oc l^^W^x,. Since the 
leptonic tensor, l^y, is known, the cross section provides us with structure information about the 
target nucleon through the hadronic tensor, W^v> which is decomposed into symmetric unpolarized 
and antisymmetric polarized parts: 

W^'^Hx,Q') = (-^^V ^) F,(x,e^) + {P^ - ^^q^^ {P^ - ^,^) (1.3) 
=/£'^^''% {^-^[g,{x,Q')+g,{x,Q^))-'LE^g,{x,Q')^ , (1.4) 
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with kinematic variables defined as V = q-P, = — M^, andx = Q'^/2v, and = \q'^\. The unpo- 
larized structure functions are Fi (x, Q^) and F2{x, Q^), and the polarized, g\ {x, Q^) and g2{x, Q^). 
Their moments are described in terms of Wilson's operator product expansion: 

2 f'dxx''-'F,{x,Q') = £ c'fl{x"),{^) + 0{\/Q^), 



dxx"-^F2{x,Q') = £ 4^(x"),(Ai) + 0(l/Q'), 
2 dxx''g,{x,Q') = £ eS^](x")A,(At) + 0(l/G2) 

Jo 



q=u,d 

dxx"g2ix,Q') = £ [4„J,?(m)-24„(x")a,(m)] +0(1/2^), (1.5) 

q=u,d 

where we suppressed the {}X^ /Q^ ,g{}i)) dependence of the perturbatively known Wilson coeffi- 
cients, c\, C2, ei, and e2- The moments, {x")q{iJ.), {x")Aq{lJ.) and d„{n) are calculable on the lattice 
as forward nucleon matrix elements of certain local operators. 

In addition, the tensor charge, (1)5^(/^), which probes transverse spin structure of nucleon, 
is beginning to be reported by experiments [9, 10]. This quantity is also calculable on the lattice 
much like the way the DIS structure function moments are calculated. 

In the following we report our lattice numerical calculations of the four iso- vector form factors, 
Gv{q^), Griq^), GA{q^), and Gp{q^), and four iso-vector moments of the structure functions, 
the quark momentum fraction, (x)„_^/(/x), the helicity fraction, (x)Au-Af/(M)> the tensor charge, 
(l)5«-5t/(At)' twist-3 coefficient d\ of the g2 polarized structure function. All but d\ are fully 
non-perturbatively renormalized. 



2. Formulation 

We use the standard proton operator, B = eabc{u^Cysdb)uc to create and annihilate proton 
states. We Gaussian-smear this operator for better overlap with the ground state with both zero and 
finite momentum. A Gaussian radius of 7 is chosen after a series of pilot calculations. Since the 
up and down quark mass are degenerate in these calculations, the isospin symmetry is exact. Thus 
we often take advantage of rotation in isospin space through Wigner-Eckhart theorem from proton- 
neutron matrix elements, (neutron|0|proton), to proton-proton ones, (proton|C?|proton). This is 
of course a well-known good approximation. We project the positive-parity ground state, so our 
two-point proton function takes the form 

C2pt(0 = £f^^) (B/3(?sink)Ba(?source)), (2.1) 
a,i3 V ^ / aj3 

with t = fsink — ^source- We insert an appropriate observable operator 0{q, t) at time T, fsource < T < 
fsink, and possibly finite momentum transfer q, to obtain a form factor or structure function moment 
three-point function, 

Cl^{t,T.q) = Y.r„p{Bp{t,inMq,t)Ba{tsou,ce)), (2.2) 
a,j5 
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rrifa #ofconfig.'s meas. interval 



sources 



(GeV) MN (GeV) 



0.005 932 10 

0.01 356 10 

0.02 98 20 

0.03 106 20 



4 
4 
4 



4 



0.33 1.15 

0.42 1.22 

0.56 1.39 

0.67 1.55 



Table 1: Number of gauge configurations and pion and nucleon mass. 



with appropriate projection, F 



l+7r 



, for the spin-unpolarized, and F 



l + 7r 



iysYk^k^A, for 



2 



the polarized. Ratios of these two- and three-point functions give plateaux for < r < t that give 



we need to take care of extra kinematics. We renormalize structure function moments by Rome- 
Southampton RI-MOM non-perturbative renomalization prescription [6]. The good continuum-like 
flavor and chiral symmetries of domain-wall fermions are very useful here in eliminating unwanted 
lattice-artifact mixings that are present in many other fermion schemes. For form factors this 
provides natural lattice renormalization through the vector charge. See our earlier publications, ref. 
[7], and references cited there in for further details. 

The RBC-UKQCD joint (2-i-l)-flavor dynamical DWF coarse ensembles [1] are used for the 
calculations. These ensembles are generated with Iwasaki action at the coupling j3 =2.13 which 
corresponds to the lattice cut off of about c?^' = 1.73(3) GeV. There are two lattice volumes, 
16^ X 32 and 24^ x 64, corresponding to linear spatial extent of about 1 .8 and 2.7 fm. The dynamical 
strange and up and down quarks are described by DWF actions with the fifth-dimensional mass of 
Mj = 1.8. The strange mass is set at 0.04 in lattice unit and turned out to be about twelve percent 
heavier than physical including the additive correction of the residual mass, wiies = 0.003. The 
degenerate up and down mass is varied at 0.03, 0.02, 0.01 and 0.005, con^epsonding to pion mass 
of about 0.67, 0.56, 0.42 and 0.33 GeV and nucleon mass are about 1.55, 1.39, 1.22 and 1.15 GeV. 
We summarize the gauge configurations and pion and nucleon mass of each ensemble in Table 1 . 

There are two important sources of systematic error: finite spatial size of the lattice and excited 
states contamination: Chiral-perturbation-inspired analysis of the former for meson observables 
suggests the dimensionless product, nijiL, of the calculated pion mass, niji and lattice linear spatial 
extent, L, should be set greater than 4 to drive the finite-volume connection below one percent, 
and the available lattice calculations seem to support this. While our present par ameters satisfy this 
condition, it should be emphasized that this criterion is not known sufficient for baryon observables. 
It is important to check this through the present calculations, and it is indeed an important purpose. 

One should adjust the time separation, t, between the nucleon source and sink appropriately 
so the resultant nucleon observables are free of contamination from excited states. The separation 
has to be made longer as we set the quark masses lighter. Here our previous study with two 
dynamical flavors of DWF quarks [7] help as the lattice cutoff is similar at 1.7 GeV. In Figure 1 we 
compare a nucleon observable, isovector quark momentum fraction, (x)„_rf, at the up/down mass 
of 0.02 where pion mass is about 500 MeV and nucleon 1.3 GeV. A cleai^ systematic difference is 
seen between the shorter time separation of 10, or about 1.16 fm, and longer 12, or 1.39 fm: the 



the bare lattice matrix elements of desired observables: e.g. {O) 



bare 



C2pt(0 



at q'^ = 0. At finite 
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Figure 1: A nucleon observable, isovector quark momentum fraction, from RBC 2-flavor dynamical 

DWF ensemble with m„^fl ~ 0.02 [7], with source-sink separation of 10 and 12: a clear systematic difference 
is seen. The shorter source-sink separation is not manifestly free of excited-state contamination. 



former averages about 0.24 while the latter about 0.20. Although the former appears statistically 
better, as the two-point function that provides the normalization has not decayed much yet, it is not 
manifestly free of excited-state contamination that is expected more at shorter time separation. The 
latter suffers statistically, as the two-point function that provides the normalization has decayed, 
but is certainly freer of excited-state contamination than the former: Of the two, we have to choose 
the latter. If larger statistics is required, then we have to collect as much statistics as necessary for 
the latter. Ideally we should collect more statistics at even larger source-sink separation to confirm 
the excited-state contamination has been eliminated by the separation of 12. 

Since in the present ensembles we explore much lighter up/down quark mass than this two- 
flavor example, with pion mass as light as 330 MeV and nucleon 1.15 GeV, the source/sink time 
separation must be chosen longer: longer than 1 .2 fm at least. In Figure 2 we present the nucleon 
effective mass at the lightest up/down quark mass, nifa = 0.005, in the present (2-(-l)-flavor dynam- 
ical DWF calculation. The nucleon signal begins to decay ait = 12, or about 1.4 fm: this is about 
longest distance we can choose without losing the signal, and hence about as free of excited-state 
contamination as we can go. As is shown in Figure 3 the bare three-point function signals for quark 
momentum and heUcity fractions for this source-sink separation of ? = 12 are acceptable. Thus we 
choose the source-sink separation of 12, or about 1.4 fm for the present work, as the two-flavor 
example suggests strongly that anything shorter would suffer from excited-state contamination. To 
confirm this choice is sufficiently long to eliminate the excited-state contamination, we need to 
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Figure 2: Nucleon effective mass at quark mass nifa = 0.005 in the present calculation with (2-i-l)-flavor 
dynamical DWF quarks. The nucleon signal begins to decay at time separation 12. By which time we are as 
free of excited state as possible while maintaining the nucleon ground state. 
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Figure 3: Bare three-point function for quark momentum fraction, (x)„_j, (left) and helicity fraction, 
{x)ai,-m, (right) for mja = 0.005 and 0.01. 

collect more statistics with a longer source-sink separation. 

3. Axial charge 

This part has been published [8]: here we repeat the essence for the readers' convenience. 
The isovector axial charge is the value of the axial form factor at zero four-momentum transfer, 
gA = Ga{0). Since the DWF action preserves the chiral symmetry, the quark vector and axial 
currents share a common renormalization up to higher-order deiscretization errors, 0((mya)^). 
Hence the lattice ratio, (gA/,?v)'^""^'^> to the vector charge, is akeady renormalized and is directly 
comparable with the experiment. In the left pane of Figure 3 we summarize the results for this 
ratio from the current RBC/UKQCD (2-i-l)-flavor dynamical DWF ensembles on large (2.7 fm) 
and small (1.8 fm) volumes, all plotted against the calculated pion mass squared, m|. 
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Figure 4: Left pane: Lattice ratio, sa/sv, plotted against pion mass squared, for the current 
RBC/UKQCD (2H-l)-flavor ensembles, large (24^ = (2.7fm)^ circles) and small (16^ = (L8fm)^ squares) 
lattice, and old RBC 2-flavor (diamonds) ensembles. 
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Figure 5: Left: Linear fit to the three heaviest points from the larger, 2.7 fm, volume, of calculated axial 
charge, sa/sv- Right: an example of simultaneous fit with finite-size effect to all the results from both large 
and small volume results. 

In the larger, 2.7 fm, volume of the present calculation, the three heavier quark mass points at 
mfa = 0.01, 0.02 and 0.03 do not show much dependence on quark mass and are close to the exper- 
imental value. However at the lightest quark mass, nifa = 0.005, the value deviates significantly 
downward and away from the experiment. In the smaller, 1.8 fm, volume, the central values all 
deviate downward from the corresponding large-volume results, albeit with large statistical errors. 
This suggests the downward shift is caused by the small finite spatial volumes of the lattice. It is 
likely such a finite-size effect is enhanced at lighter quark mass. 

The same ratio in the old RBC 2-flavor dynamical DWF calculation at a similar cutoff and 
with a lattice size comparable to the current smaller volume, 1.9 fm across, show similar trend that 
the heavier two point stay close with each other and with the experiment, and the lightest point 
deviating downward away from the experiment. These small- volume results suggest the downward 
shift away from the experiment may be caused by the finite-size effect. Indeed if we plot these 
values against another variable, the dimensionless product, m^^L, of the calculated pion mass, niji, 
and the lattice linear spatial extent, L, then the points align on a monotonically increasing curve, 
as is shown in the left pane of Figure 3. In other words, the quantity scales with the variable mj^L. 
Thus we conclude the downward shifts are finite-size effect. 

We also note that an earlier quenched RBC calculation [11] saw similar finite-size effect. 
However in these quenched calculations the effect was not detected till the linear spatial extent 
of the lattice was made shorter than 2.4 fm. It appears that the finite-size effect is enhanced by 
dynamical quai^ks. 

A similar scaling in this variable, m^iL, is obtained [8] if we replot 2-flavor Wilson fermion 
calculations of the same quantity [12, 13, 14], corroborating our discovery of the m^rL-scaling here. 



In light of this unexpectedly large finite-size effect, apparently due to light dynamical quarks, 
we need much larger volumes to accurately calculate the axial chai^ge, at least, and likely other 
elastic form factors. Since the finite-size effect grows larger than 1 % at riijiL 5, we would need 
a 5 fm across lattice for moderately light 200 MeV pions and 7 fm for realistic pions. 
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On the other hand we can exploit the finite-size effect in extrapolating to the physical pion 
mass. The simplest way is to drop all the results where finite-size effect is seen and use only the 
three heaviest points in the large volume which comfortably accommodate a linear fit: we obtain 
a value 1.17(6) for gA/gv at physical pion mass, as is summarized in the left pane of Figure 5. 
Alternatively we can fit with an additional finite-size effect term such as exp(— m;i:L) to the linear 
form and obtain values such as 1.20(6) as is described in the right pane of Figure 5. Detailed form 
of such finite-size effect term does not matter at the present level of calculations as our dynamic 
range in the scaling variable, nijcL, is limited, yet the variance among them yield an estimate of the 
associated systematics: thus we quote a value gA/gv = 1 -20(6)5131. (4)syst.- 



4. Momentum dependence of form factors 



Let us first discuss the vector current quantities, the vector and induced tensor form factors. 
The vector charge, gv = Gv{0), provides a check for the vector and axial current renormaliza- 
tion, Zy = Za = l/gv- Since the vector current is conserved, the vector charge is not affected 
by excited-state contamination, and hence is calculated very accurately. The obtained value of 
Zy = 1/1.3918(19) =0.7185(10) in the chiral limit agrees very well with an independent calcula- 
tion, Za = 0.7161(1), in the meson sector [1] upto small 0{mfaY discretization error. 

With the vector charge providing the normalization, we find the vector form factor, Gv{q^), 
is fit well by the conventional dipole form, gy /{I + {r1)q^ / 12). The single fitting parameter, the 
dipole coefficient, (r,), is the isovector Dirac mean-squared radius. The induced tensor form 
factor, Griq^), is also fit well by the dipole form, F2{0) / {Imp^) / {\ + {r\)q^ /\2), (see Figure 6.) 
Here the strength, ^2(0), is a fitting pai^ameter, and provides the isovector anomalous magnetic 
moment of nucleon, /Xp — /i„ — 1 . It agrees reasonably well with the experiment. The second fitting 
parameter, the dipole coefficient, {r^), is the isovector Pauli mean-squared radius. 
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Figure 6: Left: Dipole fit for the induced tensor form factor. Right: Resultant anomalous magnetic moment, 
Hp — H,, — I, of nucleon. The linear fit exclude the lightest point (striped circle). It extrapolates to within 
reasonable neighborhood of the experiment. 
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Figure 7: Dirac and Pauli root mean-squared radii of nucleon. Linear fits exclude the lightest points (striped 
circles). Though they increase linearly toward the lighter physical pion mass, the expected divergent chiral 
behavior is yet to be seen. Whether the lowest order heavy-baryon chiral perturbation model can describe 
such divergence is not clear in the Dirac radius, while it may still be consistent with the Pauh radius. 



The dependence of the isovector Dirac, (rj), and Pauli, (r^), mean-squared radii on the 
up/down quark mass is summarized in Figure 7. They both increase toward decreasing pion mass 
and are easily fit by linear forms, yielding extrapolations to physical pion mass that undershoot 
experiments. We do not see any sign of chirally expected divergent behavior in the chiral limit 
yet. Whether any of the chiral-perturbation-inspired models can describe these mean-squared radii 
caculated on the lattice remains to be seen: the lowest-order "heavy-baryon chiral perturbation" 
seems inconsistent with our Dirac radius result but may be consistent with the Pauli one. 

Let us now turn our attention to the axial-current form factors: The axial vector form factor 
admits single -parameter dipole fit (see the left pane of Figure 8). The dipole coefficient, (r^), 
stays constant for large values of the scaling variable m^^L about 40 % smaller than the experiment 
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Figure 8: Dipole fits to axialvector form factor (left) and resultant root mean-squared radius in comparison 
with other lattice calculations (right). Green dashed lines indicate experiments. The calculated mean-squared 
radius is well below the experiment but shows scaling in nij^L similar to the axial charge. 
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Figure 9: Left: Pion-nucleon coupling, gjiNN, extracted from near the pion pole, ~ —ni^ , of the calculated 
induced pseudoscalar form factor. The linear fit excludes the deviating lightest point, and yet undershoots the 
experiment. Right: muon-capture coupling gp = m^Gp{q^ — 0.88m^), also extracted from the calculated 
induced pseudoscalar form factor The linear fit again exclude the slightly deviating lightest point which 
would not affect the agreement with the experiment too much. 



and independent of the number of dynamical flavors, and then deviates away further downward at 
mj[L 5 and smaller where only dynamical results are available. This behavior corroborates the 
axial currents suffer from large finite-size effect. Here at finite momentum we do not know yet if 
the effect is enhanced with dynamical quarks as there is no quenched calculation with light valence 
quarks with nij^L < 6. 

From the induced pseudoscalar form factor, Gp{q^), we can extract the venerated pion-nucleon 
coupling, gnNN, from the residue at the pion pole, Gp{q^) ~ 1F„gjiMN / {'f' + fn\), ~ —tn\. As is 
seen in the left pane of Figure 9, the results for the three heavier up/down quark mass linearly ex- 
trapolate well within two standard deviation of the experimental value. The lightest point, however, 
shows downward shift away from the linear fit and the experiment, much like the axial charge. In- 
deed this coupling, g^NN, and the axial charge, gA, are related by the Goldberger-Treiman relation 

[15], q'Gp{q')-2m^GA{q^) = S^f^^f- . 

Another interesting experiment on the induced pseudoscalar form factor is muon capture. We 
extract from our calculation the effective coupling, gp = m^Gp{q^ = 0.88m^), for this process. 
As is shwon in the right pane of Figure 9, our results compare favorably with the recent MuCap 
result [16]. Here again we see the lightest point deviating downward, away from the linear fit to 
the three heavier points and the experiment. But the effect is less significant and does not affect the 
seemingly good agreement with the experiment. 



5. Moments of structure functions 



Let us first discuss the lattice ratio, (x) u-d / {x)au-m, of the isovector quark momentum fraction 
to the helicity fraction. Much like the form factor ratio, gA/gv, the momentum fraction, (x)„_(/, 
the first moment of the F\ 2 unpolarized structure functions, and helicity fraction, (x)ak-Ac/> the 
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Figure 10: Naturally renormalized ratio, {x)u-ci/ {x) Au-Ad, of the isovector quark momentum fraction to the 
helicity fraction. 



first moment of the polarized structure function, share a common renormalization as they are 
related by a chiral rotation and the DWF preserve the chiral symmetry. Thus the ratio calculated 
on the lattice is directly comparable with the experiment. The results of our calculation are shown 
in Figure 10. They do not show any discernible dependence on the up/down quark mass, and are in 
excellent agreement with the experiment. This is in contrast to the similarly naturally renormalized 
ratio of gA/gv of elastic form factors which at the lightest point deviates significantly from heavier 
mass results and the experiment. This suggests the moments of inelastic structure functions such 
as the momentum fraction, {x)u-d, and helicity fraction, {x)Au-Ad, may not suffer so severely from 
the finite-size effect that plagues elastic form factor calculations. To clarify this, we are performing 
follow up calculations at the smaller volume of 1.8 fm. 

Next we discuss the absolute values of the isovector quark momentum fraction, {x)u-d- This is 
the first moment of the unpolarized structure functions, Fi and F2. It is now fully non-perturbatively 
renormalized and run to MS at 2 GeV, Z^^(2GeV) = 1.15(4). The results are shown in the left 
pane of Figure 1 1. On the larger, 2.7 fm, volume, the three heavier points stay roughly the constant 
which is about 70 % higher at ^ 0.26 than the experiment, about 0.15. This behavior is not so 
different from old RBC quenched results [17] wtih similar up/down quark mass. The lightest point 
shows a sign of deviation away from this heavy constant behavior, but in contrast to the form 
factor deviations, is trending toward the experiment. Since lighter quark can more easily share its 
momentum with other degrees of freedom, this trending toward the experiment may well be a real 
physical effect: It is not necessarily be affected by the finite spatial size of the lattice. To clarify, 
we are performing follow up calculations at mja = 0.01 on the smaller volume of 1.8 fm, with mj^L 
smaller than the nifa = 0.005 on 2.7 fm volume. If this follow up calculation stays with the higher 
value, ~ 0.26, then this moment is not so seriously affected by the finite-size effect. 

In the Figure we also compare our results with the LHPC/MILC mixed action calculations 
[18]. Their values are significantly lower, by about 20 %, than ours. The difference is likely either 
caused by their use of only perturbative renormalization or the significantly shorter source/sink 
separation in time, or both. Another possibility is the non-unitary matching of the DWF and rooted- 
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Figure 11: Absolute value of the quark momentum fraction, (x)„_^/, fully non-perturbatively renormaUzed 
(left), and the helicity fraction, (jc)AH-Arf^ also fullly non-perturbatively renormalized (right). 



staggered fermions. 

The isovector quark helicity fraction, {x)^u-Ad, is also fully non-perturbatively renormalized 
and run to MS at 2 GeV. The renormalization agrees very well with the momentum fraction one, 
2MS^2GeV) = 1.15(3). As is shown in the right pane of Figure 11, the observable exhibits very 
similar behavior to the momentum fraction, as can be expected from the near constancy of their 
ratio: the three heavier points stay roughly the constant which is about 70 % higher than the ex- 
periment. This is again not so different from old RBC quenched results [17] wtih similar up/down 
quark mass. The lightest point then again shows a sign of deviation away from this heavy constant 
behavior, and again in contrast to the form factor deviations, is trending toward the experiment. 
Again this trending toward the experiment may well be a real physical effect, or it can be affected 
by the finite spatial size of the lattice. The situation may be clarified by the on-going small- volume 
calculations. 

Again the LHPC/MILC mixed action results are significantly lower than ours. And again the 
difference is likely either caused by their use of only perturbative renormalization or the signifi- 
cantly shorter source/sink separation in time, or their combination. Yet again another possibility is 
to the non-unitary matching of the DWF and rooted-staggered fermions. 

The results for the isovector tensor charge, {l)sii-Sd^ are presented in the left pane of Figure 
12. These are fully non-perturbatively renormalized and run to MS at 2 GeV, with Z'^^(2GeV) = 
0.783(3). These can provide a prediction since the experiments are yet to report a number. If we 
fit the heavy three points with a constant we obtain a value of about 1.10(7). Alternatively if we 
linearly extrapolate the two lightest points we would obtain about 0.7. Since we do not understand 
the finite-size coiTcction nor chiral extrapolation, these values are about as good as we can make. 

In the right pane of Figure 12 we present the twist-3 moment, di, of the polarized structure 
function gi- Unlike the other structure functions we discussed in the above, this one is not renor- 
malized yet. Our interest here is in whether the perturbatively obtained Wandzura-Wilczek relation 
[19] holds. From the smallness of the values obtained, we conclude it holds. We note we have the 
lightest mass points deviating from the linear trends set by the heavier points. 



12 



(2+] }-flavor dynamical dwf nucleon structure 



—1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 1— 



<i> 



8u-8d 



0.05 r 
0.04 '■ 
0.03 '- 
0.02 '■ 
0.01 '■ 
'- 
-0.01 ■ 



Unrenormalized 



.lat 



• all components 
■ u components 
^ d components 



0.1 



m:[GeV 



-0.02 L 



0.1 



0.2 

m^[GeV" 



Figure 12: Absolute value of the tensor charge, {l)du-Sih fully non-perturbatively renormaUzed (left), and 
the twist-3, di, moment of the g2 structure function, not renormalized (right). 

6. Conclusions 

We calculated all the four nucleon isovector form factors, vector, induced tensor, axialvector 
and induced pseudoscalar, and some low moments of isovector nucleon structure functions on the 
lattice with (2+l)-fiavor dynamical DWF quarks. The lattice cutoff is about 1.7 GeV, and two 
spatial volumes, one 2.7 fm across and the other 1.8 fm, are investigated. The dynamical strange 
quark mass is set about 12 % heavier than physical, and the degenerate up/down quark mass is 
varied corresponding to pion mass of about 0.67, 0.56, 0.42 and 0.33 GeV and nucleon mass about 
1.55, 1.39, 1.22 and 1.15 GeV. We carefully chose the nucleon source and sink separation in time 
as long as possible, at about 1 .4 fm, in order to eliminate excited-state contamination as much as 
possible while maintaining the signal. We suspect anything shorter would result in unacceptable 
contamination for the observables. 

A big finite-size effect is seen at the light up/down quark mass in the axial charge, gA/gv, a 
naturally renormalized observable with DWF quarks that preserve the chiral symmetry. The effect 
takes this observable away from the experimental value. It scales with a single variable, nij^L, and 
sets in at about mj^L ~ 5. Similar effects are seen in axialvector and induced pseudoscalar form 
factors, but are absent in vector and induced tensor form factors: This is because the axial current 
is not conserved while the vector current is. Small volume calculations no longer make sense, at 
least for axialvector and induced pseudoscalar form factors. We plan to use auxiliary determinant to 
make larger volume calculations possible. The induced tensor form factor yields a reasonable value 
for the isovector anomalous magnetic moment of nucleon, /Xp — /x„ — 1 , when linearly extrapolated 
to physical pion mass. The Dirac and Pauli isovector mean-squared radii increase toward lighter, 
physical pion mass as well, but does not show any sign of chiral divergence. 

Though similarly naturally renormalized, the ratio, {x)u-d/ {^)m-m, of the quark momentum 
and helicity fractions, is not affected by such a finite-size effect. It does not depend on light quark 
mass either and agree well with the experiment. Thus the calculations of the moments of inelasitc 
structure functions may still make sense on small volumes. Their respective absolute values, fully 
non-perturbatively renormalized, shows interesting trending toward the respective experimental 
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values at the lightest light quark mass. This may well be real physical effect driven by lighter 
quarks, in contrast to the finite-size effect in the axial-current form factors. Our on-going follow-up 
calculations at m/a = 0.01 on the smaller volume of 1 . 8 fm, with nij^L smaller than the mfa = 0.005 
on 2.7 fm volume, has a potential to resolve this question. 

In addition, the fully non-perturbatively renormalized tensor charge, {i)5u-8d been ob- 
tained. We need to understand the finite-size corrections or chiral extrapolations before transform- 
ing the values into a prediction for the upcoming experiments. 

The twist-3, di, moment of the g2 structure function is also obtained. Though yet to be renor- 
malized, its smallness suggests the Wandzura-Wilczek relation holds. 
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